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ABSTRACT 

Recent observations have discovered the presence of a Box/Peanut or X-shape structure in 
the Galactic bulge. Such Box/Peanut structures are common in external disc galaxies, and 
are well-known in N-body simulations where they form following the buckling instability of 
a bar. From studies of analytical potentials and N-body models it has been claimed in the 
past that Box/Peanut bulges are supported by “bananas”, or vivi orbits. We present here a set 
of N-body models where instead the peanut bulge is mainly supported by brezel-like orbits, 
allowing strong peanuts to form with short extent relative to the bar length. This shows that 
stars in the X-shape do not necessarily stream along banana orbits which follow the arms of 
the X-shape. The brezel orbits are also found to be the main orbital component supporting 
the peanut shape in our recent Made-to-Measure dynamical models of the Galactic bulge. We 
also show that in these models the fraction of stellar orbits that contribute to the X-structure 
account for 40 — 45% of the stellar mass. 
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1 INTRODUCTION 

There has been much recent interest in the dynamical structure of 
the Galactic bulge, following the discovery of the bimodal distribu¬ 
tion of red clump giants (RCGs) magnitudes in upper bulge fields 
(the split red clump; [Me William & Zoccali|2010[|Nataf et al.|2010[ 
[Ness et al.|2012|>. Ana lyzing 50 million RCGs from the VVV sur- 
vey, I Wegg & Gerh^dj ( |20 1 3^ were able to measure directly the 3D 
density distribution of the dominant bulge stellar population, find¬ 
ing a very pronounced Box/Peanut bulge (B/P bulge) structure ex¬ 
tending about 2kpc along its major axis. B/P bulges are found in 
about half of the external edge-on disc galaxies ( jLlitticke, Dettmarj 
|& Pohlen|2000j) and are the f ocus of several recent and on-going 

studies ( [Williams et aH201l||Fabricius et al.|20T^ [Watcher et al.| 

|2014[[^del et al.|2015| ). B/P bulges are common in N-body sim¬ 
ulations of barred discs where instability of the bar leads to one 
or more buckling events, creating a prominent peanut shape in the 
inner part of the bar ([Combes & Sanders |1 98 l[|Raha et al.|1991[ 

[Martinez- Valpuesta, Shlosman & Heller|2006^. 

The orbital structure of bars was first studied in two dimen¬ 
sions (for a review see [Contopoulos & Grosb0l||198^ , showing 
that bars are largely made out of regular orbit families trapped 
around stable periodic orbits. Later, studies of the 3D orbital struc¬ 
ture were carried out, focusing mostly on the stable periodic orbit 
families. The third dimension introduces instability in the dynamics 
and leads to bifurcations of stable periodic planar orbits. Pioneer- 
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ing work by [Pfenniger & Friedli| ( [199f[ found that the main orbit 
family of planar bars (the x\ orbits) becomes vertically unstable 
in certain regions and bifurcates to several 3D orbital families in¬ 
cluding the xivi family, xivi orbits are also called banana orbits 
because of their banana shape when seen side-on. These banana or¬ 
bits have since been considered as the backbone of B/P bulges (eg. 
Pfenniger & Friedlij 1991 [[Martinez- Valpuesta, Shlosman & Heller[ 

2006^ , even though other higher order resonant orbit families could 
also build a peanut shape bulge ( [Patsis et al.|2002| ). 


In this letter we show that the banana orbits are not necessarily 
the main building block of B/P bulges, challenging the picture of 
stars streaming along the arms of the X-shape. We describe a new 
family of brezel-like orbits with which a peanut bulge can form 
at smaller radii relative to the bar than is possible with banana or¬ 
bits. We find that this orbit family dominates the peanut shape of 
the Galactic bulge in both our N-body models (Section [^ and our 
Made-to-Measure models (Section [^ for the Galactic bulge from 
[Portail et~n?] ( [2015| . We also present an orbit-based characteriza¬ 
tion of the X-shape structure of the Galactic bulge and show that in 
our models, the Galactic X-shape structure accounts for 40 — 45% 
of the stellar mass of the bulge. 
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2 THE ORBITAL STRUCTURE OE B/P BULGES IN 
N-BODY MODELS 

2.1 N-body models of B/P bulges 

We analyze N-body models of barred discs evolved from a near 
equilibrium stellar disc embedded in different live dark matter ha¬ 
los. During this evolution the disc naturally forms a bar which 
rapidly buckles out of the Galactic plane and creates a B/P bulge 
( jCombes & Sanders|1981[|Raha et al.|1991^ . In this work we use 
the initial N-body models M80, M85 and M90 already presented in 
[Portail et aL] ( [2015^ . These three models differ by their dark matter 
fraction in the bulge, ranging from 30% for model M80 to 12% for 
model M90. Using the definition of |Debattista & Sellwood|p000t 
where M is the ratio between corotation radius and half-length of 
the bar, they span the complete range of reasonable values for the ^ 
parameter ( |Elmegreen|1996^ , from a slow bar for M80 = 1.8) 

to a fast bar for M90 =\ .08). 

We assume a distance to the Galactic centre of Rq = 8.3kpc 
( [Reid et al.|2014[[Qiatzopoulos et al.|20T5t and place the bar at an 
angle of 27° with respect to the line-of-sight towards the Galactic 
center ( [Wegg & Gerhard|2013[ ). We scale the models to the Milky 
Way by placing the end of their long bar component at I = 27° as 
seen from the Sun’s location ( [Hammersley et al.|[2000[ [Cabrera- [ 
[Lavers et al.|2007^ , resulting in a bar half length of 4.6kpc. In the 
following we use the (x,y,z) coordinates in the rotating frame of 
the bar where x is along the major axis of the bar and z is the axis 
orthogonal to the galactic plane. 


2.2 Orbit classification 


Most studies of the orbital structure of B/P bulges focus on study¬ 
ing families of periodic orbits in rotating potentials, either analyti- 
cal (e.g. [Skokos, Patsis & Athanassoula|2002[) or N-body ([Pfen- 


[niger & Friedli[ r991[ but see also [Harsoula & Kalapotharakos 

( [2009^ ). In this paper, we study the orbits of the particles in the 
N-body potential frozen well after the formation and buckling of 
the bar, which is then assumed to rotate at a constant pattern 
speed. We integrate the orbits using a drift-kick-drift adaptive leap¬ 
frog algorithm and use frequency analysis to classify those that 
build the B/P bulge. We consider all stellar particles in a box of 
±4kpc X ±1.5kpc X ±1.5kpc in the bar frame (~ 7 x 10^ parti¬ 
cles) and integrate the orbits for 50 dynamical times, where the dy¬ 
namical time is defined as the time necessary to complete a circular 
orbit at 2kpc. We record 100 positions of the selected particles per 
dynamical time and construct this way a time series for the particle 
coordinates. For each particle we compute the Fast Fourier Trans¬ 
form of its coordinate time series and identify the main frequency 
of each coordinate as the frequency corresponding to the highest 
spectral peak. We do this for the Cartesian coordinates x and z in 
the bar frame as well as the cylindrical radius r and consider the 
ratios of the main frequencies fr/ fx and fz/ fx- 

In our models we find two main groups of particles: bar 
particles and disc particles. Bar particles are identified by their 
frequency ratio fr/fx C 2 ±0.1. They have an elongated shape 
along the bar major or intermediate axis. Disc particles are defined 
through fr/fx 0 2 ± 0.1. They do not support the bar shape and do 
not show any prominent peanut shape. As we focus on the peanut 
shape of the bar, we exclude the disc particles from further analysis. 

In order to identify the orbital components of the peanut struc¬ 
ture we focus on the frequency ratio of the side-on coordinates, 
fz/fx- The top plot of Fig. 0 shows the distributions of the fre¬ 
quency ratio fz/fx for our three N-body models. As fz/fx is almost 


M80 M85 M90 



Figure 1. Fraction of the selected stellar particles belonging to orbital 
classes A-F defined through their frequency ratio fz/fx as shown in the 
figure. Top: for the three N-body models (see Section; bottom: for the 
models fitted to the Milky Way bulge (see Section[^. The colored dots in¬ 
dicate the fraction of mass in the orbital classes while the solid lines show 
the continuous distribution of the frequency ratio fz/fx- The banana orbits 
are located in the peak atfz/fx = 2. 


entirely contained between 1.5 and 2.0, we make a low resolution 
orbit classification by splitting the particles into 6 classes A to F, 
corresponding to equally spaced bins in fz/fx centreed on 1.5, 1.6, 
1.7, 1.8,1.9 and 2.0. The colored dots in Fig.[2indicate the fraction 
of stellar mass belonging to each orbital class in the three models. 

Fig.[^shows the side-on projections of these 6 orbital classes 
identified in model M85. All these classes display a strong peanut 
shape on the side-on projection. We find that orbits with larger fre¬ 
quency ratio fz/fx have larger radial extent. The total B/P structure 
of the model is then the sum of all these embedded peanut shapes. 
All orbits in any one class are very similar to each other and the 
white curves in Fig.[^show the side-on trajectory of a typical sam¬ 
ple orbit of each class. 


2.3 Which orbit classes dominate the X-shape? 

The orbital structure of 3D bars has been studied both in analyti- 
cal potentials by[Contopoulos & Barbanis[ ( [1994) and [Skokos, Pat-[ 
[sis & Athanassoula[p002^ , and in self gravitating N-body models 
by [Combes et al.[ ( [1990^ and [Pfenniger & Friedfij ( [1991[ ). These 
authors found that the dominant xi orbit family of 2D bars may 
become vertically unstable when introducing the third dimension. 
Instability of planar xi orbits results in the birth of many different 
3D periodic orbits families, called the xi tree. The simplest mem¬ 
ber of the xi tree is the xivi family, also called “banana” orbits by 
[Pfenniger & Friedli] ( |1991[ ) because of the shape of their side-on 
trajectories. Banana orbits have been found in many N-body or an¬ 
alytical models and have since been considered as the main build¬ 
ing block of B/P bulge ( [Combes et al.[[1990[ [Pfenniger & Friedli| 
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Figure 2. Side-on projections of the six orbital classes identified by fre¬ 
quency analysis of model M85. The white curves show the side-on trajec¬ 
tories of typical sample orbits belonging to each class. The banana orbits 
identified by |Pfenniger & Friedli| ( (1991) are in class F. 



X [kpc] 


Figure 3. Peanut shape (top) and X-shape produced by unsharp-masking 
(bottom), for model M85. The colors indicate the orbital class which con¬ 
tribute the highest surface density at each point of the side-on view. 


shape refers to the shape of the isophotes in the side-on projection, 
while X-shape refers to the typical X-feature that is revealed when 
applying unsharp-masking techniques to the side-on projection of 
the B/P feature. Determining the X-shape is a common practice to 
study external galaxies and can be performed in several ways, us¬ 
ing a median fi lter as in [Bureau et al.| ( [2006| ) or removing a fit of 
the image, as in |Li & Shen| ( |2012| )~ 

Fig. shows in black contours the B/P shape of the N-body 
model M85 and its X-shape obtained after removing the median 
filtered image. In this figure, the colors indicate which orbital class 
contributes the highest surface density at each point of the side-on 
view of the bar. The B/P feature, which is the sum of the different 
orbital classes shown in Fig. appears mostly boxy with no well 
defined length. On the contrary the X-shape is clear and contained 
inside 2kpc. Only classes A, B, C and D contribute significantly 
to the X-shape in this model. The banana orbits, the lowest energy 
member of the x\ tree, dominate the light only at radii larger than 
2.5kpc and contain only a small fraction of the stellar mass in the 
bulge as shown in Fig.[T] Therefore, in this N-body model, the pic¬ 
ture that the X-shape is the result of stars trapped around banana or¬ 
bits streaming along the arms of an X is incorrect. The X-shape in 
this model consist mostly of similar non resonant orbits whose par¬ 
ent orbits we call hrezel orbits, which build an X-shape at shorter 
radii. Brezel orbits are described in Section lT^ in the context of the 
Galactic bulge. 


1991 1 [Patsis, Skokos & Athanassoula||2002| 

[Martinez-Valpuesta,! 

Shlosman & Heller|2006l|Quillen et al 120131 

». They appear at the 


vertical 2 : 1 resonance (i.e. 2 vertical oscillations per revolution), 
causing the peak at fz/fx = 2 in Fig.[^ and are therefore the main 
contribution to our orbital class F. Higher order members of the x\ 
tree can also play a role in 3D bars and eventually contribute to the 
pea nut shape, as noted by [^okos, Patsis & Athanassoula| ^2002\ 
and |Patsis, Skokos & Athanassoula| ( |2002) . However these orbits 

arise at larger energies, which means even larger radial extent. 

Identifying the building blocks of the B/P feature from our N- 
body models is not a trivial task. B/P bulges are composite struc¬ 
tures where different orbits existing at different energies and radii 
overlap and all together create the B/P feature. In our models, all 
orbital classes shown in Fig. are at some radius the main com¬ 
ponent of the 3D part of the bar. In order to clarify this, we make 
the distinction between Box/Peanut (B/P) shape and X-shape. B/P 


3 THE ORBITAL STRUCTURE OE THE GALACTIC B/P 
BULGE 

3.1 Made-to-Measure models of the Galactic bulge 

In |Portail ( |20I5| ) we constructed a set of dynamical models 
of the Galactic bulge, with different dark matter halos. We used 
the Made-to-Measure method ( [Syer & Tremaine|1996[|De Lorenzi| 
|et al.|[2007t to adapt the particle weights of the N-body models 
M80, M85 and M90 just described, to force them to reproduce the 
stellar density and kinematic measurements of the Galactic bulge. 
We constrained the stellar density using the 3D density map of 
RCGs in the bulge from |Wegg & Gerha^ ( [2013^ , which trace the 
dominant fraction of the stellar mass. This map was derived by de- 
convolution of extinction and completeness corrected line-of-sight 
magnitude distributions from the VVV survey ( jSaito etal.|2012| ) and 
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Figure 4. Top: Side-on projection of the extrapolated 3D red clump giants 


number density map originally from |Wegg & Gerha^ ( (2013) . Bottom: Non¬ 
elliptical residual fraction of the side-on projection. 


covers abox of ±2.2kpc x ±1.4kpc x ±1.2kpc. As kinematic con¬ 
straints we used the radial velocity and dispersion measurements 
from the BRAVA survey ( [Kunder et al.||2012| >. After the Made-to- 
Measure fit, the resulting three models match the data very well 
and are equally good candidates to represent the dynamics of the 
Galactic bulge, differing from each other by their dark matter frac¬ 
tion in the bulge. 

The top plot of Fig. shows the side-on projection of the 
RCGs densit y map of|Wegg & Gerha^ ( [2013| ) extrapolated in the 
midplane by |Portail et al.| ( [2015^ . The peanut shape is strong and 
mostly contained inside 1.5 kpc from the centre along the major 
axis of the bar. It can be highlighted in a more physical way than 
usually done with unsharp masking. For this we follow [Portail et al.] 
( [2015^ and compute the elliptical component of the side-on pro¬ 
jection of the map E, by fitting a family of ellipses to E with the 
constraint that the residuals 'L — 'Le are everywhere positive. The 
non-elliptical residual fraction (E —Eg)/E is plotted on the bot¬ 
tom panel of Fig.[^ It shows a large deviation from elliptical shape 
at radii as short as 1.5 kpc, about one third of the bar length. In 
these models, the banana orbits do not achieve significant height 
at radii around one third of the bar length. The lower plot of Fig. 
[T] shows the orbital compositions of our three models after fitting 
to the Milky Way data. The comparison with the upper plot shows 
that the Made-to-Measure procedure (i.e. the data) moved most of 
the mass to the orbital classes B and especially C. These classes 
display a strong peanut at radii around a third of the bar length, as 
shown in Fig.|^ 

Thus, the relation between the X-shape and the orbital struc¬ 
ture is not straightforward. X-shapes in external galaxies are of¬ 
ten off-centered, with the two arms of the X not crossing at the 
center as shown by [Bureau et al.| ( p00^ . These authors suggested 
that off-centered X-shapes could be associated with banana orbits 
given their morphological similarity. The X-shapes of our Made-to- 
Measure models of the MW bulge are also off-centered (as shown 
in Fig. 17 of jPortail et al.|2015) , even though the banana orbits rep¬ 
resent only a small fraction of the orbital structure in these models. 



X [kpc] 


Figure 5. Brezel orbit viewed face-on (top) and side-on (bottom). 

3.2 The brezel orbit 

After fitting of the BRAVA data and RCGs density, most of the stel¬ 
lar mass of bar particles are on orbits of class C. These orbits were 
already present in the initial models but in lower proportion. This 
class is mostly made of non-resonant orbits which show morpho¬ 
logical similarities with one particular orbit family. Fig. shows 
the face-on (top) and side-on (bottom) view of such an orbit that 
we call brezel orbitQ This orbit has a mirror-symmetric counter¬ 
part with respect to the y = 0 plane. It has f^l — but needs 10 
oscillations in z and 6 oscillations in x in order to close. A possible 
origin for the brezel orbit would be the 5 : 3 vertical resonance of 
the so called ximul 2 planar orbit, in the notation of jSkokos, Patsisj 
|& Athanassoulaj P002^ (i.e. the usual 2D x\ orbit but considered 
to close after two rotations instead of one). Their face-on shapes 
are also similar to the rm21 and rm22 identified by jPatsis & Kat-j 
|sanikas| ( |2014j ) (see their Fig. 2), arising from a radial bifurcation of 
the ximul 2 . A detailed study of the ximul 2 tree is needed to study 
this further. 


3.3 Tbe Galactic X-structure 

What is the stellar mass fraction associated with the X-shape in the 
Galactic bulge? In models where there is no correspondence be¬ 
tween the X-shape and a particular orbit family, this question can¬ 
not be answered easily. In these models, the X-shape only reveals 
the excess of light over some smooth background in the side-on 
projection, but is no longer directly related to any physical substruc¬ 
ture of the bulge. Therefore, using our orbit classification, we define 
the X-structure as the physical stellar component which is built by 
all bar orbits (/^//x G 2 ± 0.1) that contribute significantly to the 
non-elliptical fraction identified in Fig. specifically, by all bar 
orbits that visit regions of the side-on view where the non-elliptical 
fraction is larger than 10%. With this definition, the peanut shape of 
Fig.[^can be seen as the sum of three components: the surrounding 
disc, made of orbits that do not follow the bar; the main bar, made 
of orbits that follow the bar but do not contribute significantly to 
Z — and the X-structure, contributing to most of the deviation 
of E from elliptical density. 

^ Note that our brezel orbit is topologically different from the “prezel 
boxlet” orbit, found at the 4:3 resonance in triaxial non-rotating potentials 
( [Miralda-Escude & Schwarzschildj 198^ [Merritt & Fridman|1996| . 
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With this definition, the fraction of the stellar mass in the X- 
structure contained in a box of ±4kpc x ±1.5kpc x ±1.5kpc, is 
44%, 43% and 42%, respectively, for the fitted models M80, M85 
and M90. Restricting this computation to the box of ±2.2kpc x 
±1.4kpc X ±1.2 kpc where [Wegg & Gerhard| ( [201^ measured the 
RCGs density, gives very similar numbers, 45%, 44% and 43%, re¬ 
spectively. This large fraction is in agreement with the lower bound 
estimated in |Portail et al.| ( pOT5] ). It is mostly due to the morphol¬ 
ogy of class B and C orbits, which provide the peanut shape as well 
as a significant part of the in-plane density. 


4 CONCLUSION 

We have analyzed the orbits of three B/P bulges in N-body mod¬ 
els of bars evolved in different dark matter halos. We find that the 
B/P bulges of these models are complex structures made of the su¬ 
perposition of several peanut shapes produced by different orbits, 
that are embedded at different radii. Contrary to what is usually 
stated in the literature, these B/P bulges are not mainly made out of 
“banana” or xivi orbits and their X-shape is not the result of stars 
streaming along these orbits as they follow the arms of the X-shape. 
Instead, in these models a strong peanut shape bulge is built from a 
family of brezel-like orbits, with a typical extent of about one third 
of the bar length. 

In our Made-to-Measure models fitted to Milky Way observ¬ 
ables, i.e., to the red clump giants 3D density map from [Weg'g 
& Gerhard] P013| ) and the brava kinematics data ( |Kunder et al. 
2012^ , we find that the brezel orbit family is the backbone of the 
orbits populating the Galactic Box/Peanut feature. Using our or¬ 
bital classification we propose a definition of the physical structure 
associated with the Galactic X-shape and we estimate that the frac¬ 
tion of stellar orbits that contribute to it account for 40 — 45% of 
the stellar mass of the bulge. 
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